Polarization-correlated photon pairs from a single ion 
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In the fluorescence light of a single atom, the probability for emission of a photon with certain 
polarization depends on the polarization of the photon emitted immediately before it. Here corre- 
lations of such kind are investigated with a single trapped calcium ion by means of second order 
correlation functions. A theoretical model is developed and fitted to the experimental data, which 
show 91% probability for the emission of polarization-correlated photon pairs within 24 ns. 

PACS numbers: 42.50.Ar, 42.50.Ct, 42.50.Ex 



I. INTRODUCTION 

One of the most relevant tools for characterizing the 
light field emitted by a quantum optical system is its in- 
tensity correlation function <?^(t), the most prominent 
example being the observation of antibunching [l|, 0] in 
the fluorescence of a single atom. Intensity correlation 
functions of the fluorescence light of a single ion excited 
by two light fields have been investigated before. In Q 
the pair correlation conditioned on the wavelength of the 
photons scattered by a Ba + ion (eight-level structure) 
reveal the transient internal dynamics of the ion, which 
is characterized by optical pumping and the excitation 
of Raman coherence. In a detailed description of these 
experiments [|[, it is shown that measuring the correla- 
tion function for only one polarization of one of the two 
light fields by adding a polarization filter projects the 
ion in a coherent superposition of its energy eigenstates. 
This may be considered a first signature of atom-photon 
entanglement 

Other investigations of the resonance fluorescence of 
coherently driven single atoms by means of (t) com- 
prise, e.g., a theoretical analysis of the effect of the quan- 
tized motion of a two-level atom in a harmonic trap 
on the second order correlation function [||, experimen- 
tal measurements of photon correlations revealing single- 
atom dynamics in a magneto-otical trap 0, Q and atomic 
transport in an optical lattice 0] , the theoretical demon- 
stration of nonclassical correlations in the radiation of 
two atoms [Ioj |. i.e. bunched and anti-bunched light is 
emitted in different spatial directions, and the study of 
intensity-intensity correlations from a three-level atom 
damped by a broadband squeezed vacuum (Tlj . 

Another application of photon-photon correlations in 
the resonance fluorescence from single atoms is an ex- 
periment where, using a self homodyning configuration, 
the second order correlation function was used to char- 
acterize the secular motion of a trapped ion [T^ . The 
experiment revealed the dynamics of both internal and 
external degrees of freedom of the ion's wave function, 
from nanosecond to millisecond timescales, in a single 
measurement run. 

Polarization correlations of subsequent photons enit- 
tcd from a bichromatically driven four-level atom have 



been predicted theoretically [13[ . In this proposal a four- 
level atom (J = 1/2 to J = 1/2) with degenerate Zeeman 
sub-states is coupled to a light field consisting of two com- 
ponents which are symmetrically detuned from the reso- 
nance frequency of the atomic transition. In contrast, in 
the experiment presented in this paper, the Zeeman de- 
generacy is lifted and the corresponding atomic transition 
(four-level system) is driven monochromatically. Fur- 
thermore, the effect of the coupling of the excited states 
to another manifold (J = 3/2) of four Zeeman sub-levels 
by a second light field is taken into account. 

In the context of quantum optical information tech- 
nologies, many entanglement protocols str ong ly rely on 
a projective measurement [13, [lj| [H, [13, EM, which is 
also the origin of the effect of antibunching. In fact, the 
degree of antibunching is a benchmark for single quantum 
emitters or ideal single photon sources that are suitable 
for quantum networking and communication. 

Beyond proving that a source is a pure single quantum 
emitter, correlation functions are used to characterize 
and develop applications for quantum communication. 
In [l9j . for example, the resonance fluorescence from a 
continuously excited single ion is split and recombincd 
on a beam splitter with a relative delay, thus creating an 
effective two-photon source. Here the measured correla- 
tion function reveals the quality of the mode matching at 
the beam splitter. Measuring and controlling the degree 
of second order correlation of a photon source is therefore 
a very important step in designing quantum optical tools 
in quantum information processing. This engineering of 
correlation functions is of special interest for quantum 
networks where single photons mediate information be- 
tween nodes of single atoms [2(| [2l| . 

In this context, we present a polarization selective mea- 
surement of the correlation function g^ (t) of fluores- 
cence photons from a single 40 Ca + ion, which is contin- 
uously laser excited on the 4 2 S!/ 2 to 4 2 P 1 / 2 transition. 
In particular, the correlation function of emitted o~ and 
(j + polarized photons [2rj| was measured conditioned on 
the previous detection of a a~ photon, and it is shown 
that the system can be tailored such that the polariza- 
tion of a photon depends strongly on the polarization of 
a previous one. 
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FIG. 1: Level scheme of 40 Ca+ for a non zero magnetic 
field B. The 2 Si/ 2 , 2 Pi/2 and 2 D 3 / 2 states split into 

eight levels according to their magnetic quantum 
number rrij and the Landc factors gj. The splitting 
shown is not to scale. 



II. MODEL 

Figure Q] shows the relevant levels of the 40 Ca + ion. 
Continuous excitation with laser light at 397 nm and 
866 nm generates continuous resonance fluorescence. We 
consider the case of observation along the direction of 
the magnetic field and write g^(r) in terms of the tran- 
sition operators for the cr-transitions from Sx/ 2 to P1/2 
(see figure Q] for the numbering of the levels) , 



a x = |1)(4| 
o 2 = |2)(3| 



&1 = I4)(l|, 
4 = |3}(2|. 



(1) 
(2) 



If we consider the general case without polarization se- 
lective detection, the second order correlation function 
reads (H 



9 {2) (r) 
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(3) 



Using the quantum regression theorem [231 ] it can be 
shown that for all initial conditions and at steady state, 
the second order correlation function is related to the 
excited state populations by 



.,(2) 



(r) 



P33(» +P44(t) 
P33(oo) + ,944(00) ' 



(4) 



where p; rn (r) with l,m € [1,8] are the matrix elements 
of the density operator representing populations and co- 
herences of the states \l) and \m) at time r after the 
detection of a first photon. Equation [4] allows to predict 
g( 2 ) by solving the optical Bloch equations. The second 
order correlation function g^ 2 ~> (r) of the 397 nm fluores- 
cence light of a single 40 Ca + ion is thus proportional to 
the population of the ion's excited state P1/2 at time r 
whereby the initial state at r = is determined by the 
first measured photon. 



We now consider polarization-selective detection of 
photons for the situation of our experiment. Figure [5] 
shows a sketch of the four levels involved in the emission 
of blue (397 nm) photons on the P1/2 to S1/2 transition. 
The detection of a o~ polarized photon projects the ion 
into state |2) = |Si/2, rn, = 1/2} (figure l2aj) . If the ion 
is continuously illuminated by linearly polarized light at 
397 and 866 nm with polarization perpendicular to the 
magnetic field, no 7r transitions are excited. The ion effec- 
tively sees a superposition of a~ and a + polarized light. 
After the emission of a a~ photon and the correspond- 
ing projection into state |2), a subsequent excitation can 
thus only occur to state |3) = \P\/2,rrij = —1/2) by ab- 
sorbing a cr~ polarized photon. From state [3} the ion 
can either decay back to state |2) under emission of a 
second a~ polarized photon (figure [2b]) . or it can decay 
to state |1) = |Si/2,TOj = —1/2) under emission of a -n 
polarized photon. Thus the next emitted photon after 
a cr~ photon can not be a + polarized. In order to emit 
a <7 + polarized photon, the ion must first be excited to 
state |4) = |Pi/2,TOj = 1/2), which can only happen out 
of state |1) by absorbing a a + polarized photon from the 
exciting beams (figure [2c]) . The probability of detecting 
a o + polarized photon after the time r conditioned on 
the previous detection of a a~ polarized photon is there- 
fore much lower than the probability of detecting a o~ 
polarized photon for the same condition. Since the corre- 
lation function is a measure for the photon-photon wait- 
ing time distribution, the g^ derived from correlating 
<7~ with a + photons should be suppressed with respect 
to the g( 2 > derived from correlating a~ with a~ photons. 
In fact both correlations are expected to exhibit almost 
ideal antibunching, with the difference that for the first 
case the dip around t~0 is expected to be much wider. 
In other words, and taking into account also the case 
where er + and er - are exchanged, there will be a much 
stronger antibunching for the correlation of photons with 
opposite a polarization than for photons with the same 
a polarization. The same considerations are applicable 
for excitation with pure 7r-polarized light. In that case 
the correlation of subsequent photons is much stronger 
for opposite than for equal a polarization. 

Analogous to equation 2] the second order correlation 
functions for a~ and a + light, conditioned on the previ- 
ous detection of a cr~ photon (p(0) = p- ln it = |2)(2|) are 
derived to be 



and 



P33(t) 
P33(oo) 



,944(00) ' 



(5) 



(6) 



respectively. In contrast to the general case (equation 2]) 
the conditioned second order correlation function is given 
by only one of the excited state populations at time r. 

In figure l3"al and l4"al the correlation functions calculated 
according to equations [5] and [5] arc plotted for weak and 
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(a) Projection (b) a polarized photon (c) <r+ polarized photon 



FIG. 2: Detection of a~ (b) and cr + (c) polarized photons, conditioned on the previous detection of a a~ (a) 
polarized photon. Excitation happens with er + -(7~-polarizcd light, i.e. linear polarization perpendicular to the 
magnetic field, (a) Detection of a <r~ polarized photon projects the ion into the state |2) = \Si/2,mj = 1/2). (b) To 
detect a subsequent o~ polarized photon, the ion has to be re-excited under the absorption of a o~ polarized 
photon, (c) To detect a a + polarized photon, the ion has to be re-excited under the absorption of a a~ polarized 
photon, then decay to state |1) under emission of a 7r polarized photon and then be excited to the state |4) by 

reabsorbing a a + polarized photon from the exciting beams. 



strong excitation parameters typical for our experimental 

( 2) 

setup. The upper curves, in blue, show g 1{t), while the 

C2) (2) 

lower ones, in red, represent g +(r). As expected g _(t) 

(2) 

and g + (t) show very different behaviors. For both weak 

(2) 

and strong excitation, g K _ (r) rises to high values with 
a steep slope. For weak excitation a maximum g( 2 ' of 
15.6 is reached at r = ±29 ns. For longer time differ- 
ences the correlation function falls monotonously until 

(2) 

it reaches one. For strong excitation g 1(t) reaches a 
maximal value of 2.8 at t = ±13 ns. For r > 13 ns the 
correlation function falls to 1 after 200 ns showing some 
coherent oscillations. 

(2) 

The two g + ( T ) functions, on the contrary, show a flat 
behavior on short time scales before they rise with a mod- 
erate slope to their maximum value, which is reached at 
approximately 400 ns for the weak excitation and 200 ns 

( 2) 

for the strong excitation. In the latter case g +(t) rises 
directly to 1 without a transient overshoot, while for weak 
excitation it reaches a value of 3.1 before it falls to 1 for 
r > 400 ns. 

( 2) 

The shape of g a _ (r) in figure [3a] is characterized by 
large correlation values and a long decay time (compared 
to the lifetime of the P1/2 state) to the asymptotic value 
and can be attributed to optical pumping into the D 3 / 2 
state [1]. The excitation on the S1/2 — P1/2 transitions 
is much stronger than the one on the D 3 / 2 — P1/2 tran- 
sitions. As a result, a large fraction of the population 
is transferred to the D 3 / 2 state after long delay times 
and the observed fluorescence is weak. This small flux 
of fluorescence, caused by a small steady-state popula- 
tion of the state |3), determines the normalization factor 
for long time intervals ,033(00). At short time delays af- 
ter the projection into state |2), however, i.e. during the 
first 30-40 ns, a large fraction of the population is excited 
to state |3) and the optical pumping to the D 3 / 2 states 



is negligible. Since the correlation function is the ratio 
of the population of state |3) at time r and that in the 
steady state, high values are reached, which then decay 
to one revealing the time scale of the optical pumping. 

(2) 

The characteristics of g a + (r) can be explained with an 
analogous argumentation. Here, excitation to state |4) 
during the first 30-40 ns is even smaller than the steady 
state population ,044(00) for long time delays. After the 
projection into |2), it takes several scattering events and 
therefore more time until the population of state |4) ex- 
ceeds ^44(00). As the inset of figure [3a| shows, for large 
(2) 

time intervals gr + (r) decays to the asymptotic value with 

(2) 

the same time constant as g 1{t). 

The correlation functions shown in figure 0a| are calcu- 
lated for equal excitation strength on the Si/ 2 — Pi/ 2 and 
D 3 / 2 — Pi/ 2 transitions. Consequently, the correlation 

(2) 

values of g i (r) are smaller and the decay to the asymp- 
totic value is much faster. The exciting fields arc strong 
enough to cause some damped oscillations at the general- 
ized Rabi frequency VIq = y ^397 1 2 + A| 97 in the corre- 
lation function. Due to the complex eight-level structure, 
the oscillations do not occur at only one generalized Rabi 
frequency, but each transition between the Zeeman sub- 
levels contributes a Fourier component depending on the 
intensity and detuning of the driving field. In figure [4a| 

this becomes noticeable by comparing the frequency of 

( 2) 

the strongly suppressed oscillations of g +{t) with the 
(2) 

ones from g '(t). Due to the Zeeman splitting, the de- 
tunings of the |1) to |4) and the |2) to |3) transitions 
with respect to the exciting laser give rise to different 
generalized Rabi frequencies. 

The most interesting feature in figure [3al and [4al is the 

behavior of the correlation function for times t close to 
(2) 

zero, g + (t) shows a flat plateau of values very close 
to zero for —15 ns < r < 15 ns in the case of weak and 
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FIG. 3: (a) Conditioned second order correlation functions g 1{t) (top, blue) and g +{t) (bottom, red) for weak 
excitation. The functions have been calculated with a eight level Bloch equation model for Rabi frequencies of 
fl 3g7 = 2tt • 9.2MHz, ft 866 = 2tt • 1.3 MHz, detunings A 397 /27r = -15 MHz, A 866 /27r = 5.8 MHz and a magnetic field 
of B = 3.5 G. (b) Ratio p(r) (equation [7]) of the integrals over time of the two conditioned correlation functions. 
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FIG. 4: (a) Conditioned second order correlation functions g 1{t) (top, blue) and <7^ + ; (t) (bottom, red) for strong 

excitation. The functions have been calculated with a eight level Bloch equation model for Rabi frequencies of 
fl 397 = 2tt • 20.2 MHz, 866 = 2?r • 20.3 MHz, detunings A 3 9 7 /27r = -15 MHz, A 866 /27r = 5.8 MHz and a magnetic 
field of B = 3.5 G. (b) Ratio p(r) (equation[7]) of the integrals over time of the two conditioned correlation functions. 



for —10 ns < r < 10 ns in the case of strong excitation. 

(2) 

In the same time interval g _ (t) rises very fast to high 
values for both excitation conditions. The ratio 



p(t) 



(7) 



of the number of a - and cr + -polarized photons emitted 
by the ion in r describes the purity of the polarization 



in that time window. The probability with which, if a 
second photon is detected within this time interval, the 
polarization of this second photon is er~, is easily cal- 
culated as . Figure [3b] and [4b] show p for weak and 
strong excitation conditions, respectively. The ratio di- 
verges for r — > (> 10 5 for t = 1 ns) in both excitation 
regimes, because p^(r) increases quadratically in time 
from t = 0, while pair) increases with t 4 . In other 
words, if a second photon is emitted within a short time 
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FIG. 5: Setup for the measurement of 
polarization-conditioned correlation functions. The 
fluorescence light is split in to two parts by collecting it 
with the two HALO lenses. Multimode fibers are used 
to couple the light to two PMTs. 
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interval after a first er~-polarizcd photon, then with very 
high probability (that can be chosen by the time inter- 
val) it will also be cr~-polarized. By exciting the ion 
with 7r-polarized light it is possible to generate equally 
strongly correlated photon pairs of oppositely cr-polarized 
photons. 



FIG. 6: Excitation spectrum of a single ion. The 397 

and 866 nm lasers are approximately vertically 
polarized and propagate under 90° to the quantization 
axis (see figure [5]). The solid line (red) is the calculated 
spectrum for the Rabi frequencies ^397 = 2ir • 9.9 MHz, 
il 866 = 2tt • 1.5 MHz, the detuning A 3g7 /2n = -15 MHz 
and a magnetic field of B = 3.5 G. The background is 89 
c.p.s. and the data points shown have Poissonian errors. 



III. EXPERIMENT 



Figure [5] shows the setup of the measurement which 
is described in detail elsewhere [24l. |25|. A single 40 Ca + 
ion is trapped in a linear Paul trap and continuously ex- 
cited by two lasers at 397 nm and 866 nm. The 397 nm 
laser is red detuned to provide Dopplcr cooling while the 
resonant 866 nm laser prevents optical pumping to the 
D3/2 state. The 397 nm fluorescence light is split into 
two parts by collecting it with two high numerical aper- 
ture lenses (HALO) 24]. Both lenses direct the collected 
photons to photo multipliers (PMT) through multimode 
fibers. In each beam path a A/4 plate and a polarizing 
beam splitter are used to select the respective polariza- 
tion. The arrival times of the signals from the PMTs 
are recorded with picosecond resolution by commercial 
counting electronics [23], and the correlation function is 
obtained by postprocessing the data. 



1. Calibration 

Before measurement of the correlation functions, an 
excitation spectrum, i.e. the rate of 397 nm fluorescence 
as a function of the detuning of 866 nm, was recorded. 
The spectrum that shows four dark resonances provides 
a calibration of the experimental parameters, which are 
then used as the starting point to fit the conditioned cor- 
relation functions. In figure [6] this excitation spectrum is 
plotted. The Rabi frequencies, dctunings and the mag- 
netic field which are deduced from a fit to the data are 
indicated in the caption. To account for experimental de- 



viations in the polarization angle of the excitation lasers 
from the ideal vertical polarization, this parameter was 
also varied in the fit. The angles of the polarizations of 
the two laser beams with the magnetic field, as used in 
the model shown in figure are a^gj = 0.46 • it for the 
blue laser and a^ee = 0-4 • 7r for the infrared laser (rather 
than tt/2 in the ideal case). This is compatible with the 
available control over the experimental settings and the 
quality of the optical components used. 

2. Conditioned correlation functions 

Figure [7] shows the data for the two measured condi- 
tioned correlation functions plotted in one graph. Data 

( 21 

for g I are presented in red (lower curve) and data for 
( 21 

g _ in blue (upper curve). The data are normalized to 
a long-time value of one and presented without back- 
ground subtraction. The solid lines are the g^ 1 functions 
obtained by a fit of the Rabi frequencies to the experi- 
mental data using the model discussed in section [TTJ The 
values of the background, laser dctunings, magnetic field 
as well as the polarization angles extracted from the fit 
to the excitation spectrum have been kept fixed. The 
two Rabi frequencies have been fitted to both curves at 
the same time and agree well with the experiment. The 
deviation of the Rabi frequencies between the fitted cor- 
relation functions and the fitted excitation spectrum lies 
within the statistical error [28| . The correlation functions 
are very similar to the ones for the case of weak excitation 
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FIG. 7: Blue (top): (g^l). Red (bottom): (g™). The 
solid lines (green) are the correlation functions expected 
from the theoretical model for the Rabi frequencies 
Q 3g7 = 2tt • 9.2 MHz, n S66 = 2tt • 1.3 MHz and the 
detuning A 8 66/27r = 5.8 MHz. The inset shows data and 
model for time scales up to 5 /zs. 




t (ns) 

FIG. 8: Zoom of FIG. [7J 



in figure [3] For large r (> 400 ns) both functions overlap 
and slowly decay to one. The characteristic behavior of 
large (small) correlation values and the slow decay to the 

asymptotic value for gjl ) , explained in section [ill 
are clearly observed in the measurement. 

Figure [5J shows a zoom into the region for small time 
differences. At t = both curves reach a value close 

(2) 

to zero. With increasing r, g _ rises with a very steep 
quadratic slope and reaches a value as large as 12, 

(2) 

whereas g + stays flat for ~ 5 ns, before it rises with 
a moderate slope to a maximum value of 2.9. The model 



agrees very well with the data, proving the good control 
over the creation of polarization correlated photon pairs 
that we obtained experimentally. 

The main difference between the model calculation of 
g^> in figure l3"al and experimental data in figure |S] is that 

A ( 21 

the predicted r -like plateau of g K I (r) is less pronounced 
in the measurement. Simulations with the model from 
section [TT] show that this is explained by small errors in 
the polarization of the exciting lasers and in the detec- 
tion setup. The polarizations of the exciting lasers have 
been fitted to the excitation spectrum (fig[BJ), and these 
results have been used in the model of the correlation 
functions. To achieve an agreement of the quality as it 
is shown in figure [7] and [SJ also deviations from the ideal 
polarizations in the detection were accounted for in the 
model. These deviations occur when the polarization is 
not filtered perfectly, and consequently the measured g^ 
function contain some coincidences that are caused by 
the respective orthogonal polarization. The theoretical 
model in figures [7] and [8] has been calculated using 2.5 % 
of wrong a + polarized photons for the initial detection 

(2) 

events in both curves. For q _ the conditioned detection 
of the second a~ polarized photon has an error of 5%, 
while for g^ 2 } the detection of the a + polarized photon 

( 2) 

has an error of 1.8%. This demonstrates that g { + , in 
particular its r 4 characteristic, is very sensitive to small 
polarization errors at times r close to zero. These er- 
rors are within the precision with which the polarization 
filtering was controlled in the experiment, given the low 
light levels and imperfections of the optics. Furthermore, 
there is a contribution due to the collection of light from 
a solid angle along the quantization axis. Each HALO 
lens collects 4% of the light emitted into the full solid 
angle [24|, which leads to the detection of a small frac- 
tion wrong a and n polarized photons. 

The sensitivity of the generation of polarization corre- 
lated photon pairs to polarization errors gets even more 
evident considering the purity p (equation [7]) • Figure [5] 
shows this ratio for the values of the model calculation 
that was fitted to the data from figure [7J and [5] In con- 
trast to the ideal case from figure [3bl p does not diverge 
for times r — > 0, but instead reaches a maximum of al- 
most 10 at 24 ns and falls then to 9.3 at 1 ns. This means, 
if a second photon is detected with our setup within 24 ns, 
it is cr~ polarized with 91 % probability. A practical fig- 
ure of merit also has to consider the absolute number of 
photons within the time interval r. After 24 ns the ion 
emits on average 0.07 er~ polarized photons into the full 
solid angle, out of which 8% are detected with our setup 
[2~i| . Increasing the time window will yield more photons, 
but decrease the polarization purity. The curve for the 
ideal case in figure [3b] reaches a value of 130 at 24 ns, 
suggesting that it is in principle possible to generate po- 
larization correlated photon pairs with more than 99% 
probability using this method. This could be achieved 
by reducing polarization errors in excitation and detec- 
tion. 
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FIG. 9: Purity p for the model calculation fitted to the 
measured data from figure [7] and [8] To avoid artifacts at 
times close to zero p is presented without background. 



The excitation conditions of the single ion also affect 
the efficiency of our photon pair source. For the strong 
excitation conditions of figure [4j for example, the ion 
emits on average 0.2 a~ polarized photons within 24 ns. 
Simultaneously the purity decreases, yielding a cr~ po- 
larized photon in only 96% of the cases (figure [4b| . For 
a good comparison of the source under weak and strong 
excitation conditions it is convenient to look at the ef- 



ficiency of both cases for equal photon numbers. Un- 
der strong excitation conditions the ion emits on average 
0.07 photons in 12 ns. If a second photon is detected 
within these 12 ns, it is a~ polarized with more than 
99%. Increasing the Rabi frequencies of the excitation 
beams thus allows to reduce the time interval in which 
the photon pairs are emitted while keeping photon num- 
ber and polarization purity constant. 

Summarizing, it was shown that the second order cor- 
relation function of the fluorescence light of a single 
ion can be engineered by polarization-sensitive detection. 
The functions for a + and a~ light conditioned on the 
previous detection of a a~ photon show a characteristic 
anti-bunching behavior that allows for heralding photons 
of a certain polarization in an otherwise randomly po- 
larized stream of photons, with possible applications in 
quantum optical information technology. A single, laser 
cooled ion generated polarization correlated photon pairs 
within a time window of 24 ns with an efficiency of 91%. 
Model calculations show that the presented method has 
the potential to reach an efficiency of more than 99% 
within the same time window. 
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